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Water absorption in the metal-organic framework (MOF) material Cu-BTC, up to a concentration 
of 3.5 H 2 O per Cu ion, is studied via density functional theory at the meta-GGA+U level. The 
stable arrangements of water molecules show chains of hydrogen-bonded water molecules and a 
tendency to form closed cages at high concentration. Water clusters are stabilized primarily by 
a combination of water-water hydrogen bonding and Cu-water oxygen interactions. Stability is 
further enhanced by van der Waals interactions, electric field enhancement of water-water bonding, 
and hydrogen bonding of water to framework oxygens. We hypothesize that the tendency to form 
such stable clusters explains the particularly strong affinity of water to Cu-BTC and related MOFs 
with exposed metal sites. 


Introduction 

The need to reduce greenhouse gas emissions has 
driven the search for materials to capture and sequester 
carbon dioxide at a low-enough cost to be economically 
viable. Metal-organic framework (MOF) materials offer 
great promise for the capture of CO 2 . mm Absorption 
and desorption from MOF materials via swings of tem¬ 
perature and pressure in principle requires less energy 
than do current aqueous capture technologies. [2] In ad¬ 
dition, MOF materials can be created with a wide vari¬ 
ety of pore sizes, shapes, connectivities, and topologies. 
These attributes, as well as the effects of cation substi¬ 
tution on metal sites, can in principle be used to tailor 
both the physical and chemical interactions of the MOF 
with gas molecules to optimize desirable properties such 
as CO 2 selectivity, capture rates, and costs of absorp¬ 
tion/desorption cycles, [l] 

A particularly interesting and well-studied example of 
a MOF material is Cu-BTC, also known as HKUST- 


FIG. 1: Cu-BTC viewed along (a) 100 and (b) 110 directions. 
Cu atoms are blue; O red; C gray; H small yellow. 
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1 (Ref. [3]). This material (Fig. [l]) consists of copper 
dimers linked by 1,3,5-benzenetricarboxylate C 6 O 9 H 3 
(BTC) units. The structure has a three-dimensional cu¬ 
bic framework with channels of alternating 13.3 A and 
11.1 A cuboctahedral pores connected by 6.4 A square 
windows along 100 -type directions, and 5.5 A tetrahe¬ 
dral side pockets connected to the 13.3 A pores via 3.7 
A triangular windows. [H E] 

As formed in atmospheric conditions, Cu-BTC con¬ 
tains a significant degree of water. This water can be 
removed via heat-treatment under near vacuum, ME] 
leaving a structure with partially “exposed” Cu ions0 
facing the large cuboctahedral pores, allowing for par¬ 
ticularly strong interactions with adsorbates at these 
Cu sites. These exposed ions lead to several poten¬ 
tial applications in addition to carbon capture, including 
catalysis, [B] hydrogen storage, [9j storage of other gases 
such as NO, m and gas separation. m 

Dry Cu-BTC has a CO 2 uptake of as much as 19.8 
weight % at atmospheric pressure, [H] and a high selec¬ 
tivity of CO 2 over N 2 . uma The situation of hydrated 
Cu-BTC is interesting. Up to about one H 2 0 per Cu 
site, there is theoretical and experimental evidence for 
a slight increase in CO 2 uptake. [131115] For water con¬ 
centrations greater than about one H 2 0 per Cu, a situa¬ 
tion which we call here “highly-hydrated Cu-BTC”, the 
CO 2 uptake is reduced. puna Hz] At high-enough con¬ 
centrations of H 2 O, Cu-BTC loses almost all of its CO 2 
capacity. p(4] The last fact likely prevents the use of Cu- 
BTC as a material for post-conrbustion carbon capture, 
as flue gases in coal-burning plants contain significant 
water vapor. [TBj In some MOFs, water can even break 
down the structure completely. m 

In spite of the importance of the highly-hydrated state 
of Cu-BTC for its performance, little is known about 
the structure of water in this state. Absorption isotherm 
measurements [14] 19] show that as much as 32 to 40 mol 
kg -1 of water can be absorbed into Cu-BTC. Assum¬ 
ing an ideal pore structure, this corresponds to 6.5 to 
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TABLE I: Total binding energies of water dimers and hex- 
amers (in eV). VASP results are compared with benchmark 
GAUSSIAN results. The GAUSSIAN results show the av¬ 
erage and (in parenthesis, with units of the least significant 
digit) root mean square deviation for three different highly- 
converged basis sets, [29] as tabulated in Ref. [30. Only the 
electronic contribution to binding energy is included. 

DFT code water dimer water hexamer 
VASP GGA -0.30 -2.92 

VASP meta-GGA -0.23 -2.13 

GAUSSIAN -0.23(1) -2.05(16) 


8.0 H 2 0 per Cu. An NMR study[20] suggests that the 
state of highly-hydrated Cu-BTC has one water molecule 
bound to each Cu in equilibrium with fluid water in the 
remaining pore volume. A recent X-ray powder diffrac¬ 
tion refinement, ED on the other hand, shows evidence for 
2.3 bound water molecules per Cu atom, in three differ¬ 
ent partially occupied binding sites, although it was not 
possible to determine the hydrogen positions for these 
water molecules. 

Previous electronic structure calculations and molecu¬ 
lar simulations have generally investigated the structure 
of H 2 0 in Cu-BTC only up to one H 2 0 per Cu ion, [221 
[25] in which case one water molecule binds to each Cu 
position, with the water oxygen atom (Oyy) closest to 
the Cu. In this work, we use density functional theory at 
the “meta-GGA + U” level to investigate the structure 
of water absorption in Cu-BTC up to 3.5 H 2 0 per Cu. 

Computational Methods 

First principles density functional theory calculations, 
as encoded in the VASP software (Ei E3), were used 
to calculate the relaxed configurations investigated here 
and their electronic structures. All calculations were per¬ 
formed for a primitive cell of Cu-BTC containing 156 
framework atoms plus any H 2 0 adsorbates. Because of 
the large cell, only a single k-point at the origin was used. 
The planewave cutoff was 500 eV for all calculations. Van 
der Waals forces were treated using the “DFT-D2” ap¬ 
proximation of Grimme, |28j and were included in all cal¬ 
culations. 

To study water absorption in Cu-BTC, it is crucial 
to get the interactions between H 2 0 molecules correct. 
We compared various exchange-correlation (XC) func¬ 
tionals to see if the correct binding energies of small 
water clusters could be obtained. The results for two 
such XC, the PBEsol|3T] generalized gradient approxi¬ 
mation (GGA), and the PBEsol plus RTPSS[32l meta- 
GGA functional (one that also uses the second deriva¬ 
tive of the charge density), are shown in Table [TJ and 
compared with benchmark calculations of water bind¬ 
ing energies @2] that were performed using the DFT code 


GAUSSIAN[27] with highly-converged basis sets.[25] The 
results show 30 % to 40 % water-water overbonding for 
the case of the GGA calculations, but excellent agree¬ 
ment with the GAUSSIAN results for the meta-GGA ap¬ 
proximation. The meta-GGA approximation was there¬ 
fore used for the rest of this work, despite being several 
times more computationally expensive than the GGA ap¬ 
proximation. 

As noted in previous DFT works, improvement in the 
agreement between DFT results and experiments for sys¬ 
tems containing magnetic ions can generally be achieved 
if on-site Coulomb terms are included. [33] in what is 
commonly termed the “GGA + U” approach. To de¬ 
termine the onsite Coulomb parameters to use for Cu- 
BTC, we began with the experimental structure [Ml 155] 
of Cu 2 C 03 ( 0 H ) 2 (malachite). Malachite has numerous 
structural features analogous to those in Cu-BTC, in¬ 
cluding the same species (Cu,0,C,H), Cu in the Cu 2+ 
valence state, Cu coordinated to a square of oxygen, and 
OH units which are analogous to the H 2 0 admolecules 
in hydrated Cu-BTC, but malachite is simpler to inves¬ 
tigate because of its smaller unit cell. 

Our aim was to adjust the U parameter for both Cu 
and O so as to fit the bandgap and experimental crys¬ 
tallographic structure of malachite as well as possible, 
and, assuming transferability, to use the same values to 
study Cu-BTC. We were unable to find any value of the 
malachite bandgap in the literature. Our calculations, 
however, suggested that a bandgap of 1.75 eV leads to a 
minimum of absorption at around 2.3 eV, at the charac¬ 
teristic green color of malachite, and we thus fit to this 
value. We used the malachite coordinates given by Zigan 
et al. (Ref. [34] cited in Ref. [36]) , because hydrogen po¬ 
sitions were given. To quantify the structural agreement 
with experiment, we fixed ions at their experimental posi¬ 
tions and calculated the root mean square residual forces 
on the ions. For a calculated bandgap of 1.75 eV, the 
residual forces were minimized for Cu U = 3.08 eV and 
oxygen U = 7.05 eV. The magnitude of the U value for 
Cu largely controls the calculated bandgap, by (control¬ 
ling the splitting in the Cu d levels), while including a 
nonzero U for oxygen greatly reduces the residual forces. 
Magnetism of Cu 2+ ions with their d 9 electronic config¬ 
urations was treated using spin-polarized DFT calcula¬ 
tions. The antiferromagnetic arrangement of the 2 Cu 
on the Cu-Cu dimer was found to be lower in energy 
than the ferromagnetic one, in agreement with previous 

studies. ED 

We now hydrate the Cu-BTC structure. We begin 
with the positions of the water molecule oxygens (Ow) 
determined by Wong-Ng et aZ. [21] (Fig. [2j Table |III| ) , 
all located inside the large 13 A pores. In this study, 
the Ow(I) sites (green in Fig. [2]) were found to have 
an occupancy of 0.80(2), the Oq/(II) sites (blue) oc¬ 
cupancy 0.32(1) and the Ow(HI) site (red) occupancy 
0.57(2). The partial occupancy values imply an average 
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FIG. 2: (a) Large spheres show positions of water oxygens 
(Ow) inside the large pore of Cu-BTC, according to the struc¬ 
ture refinement of Wong-Ng et al. m Sites are color-coded 
according to symmetry (green: Ow(I); blue: Ow(II); red: 
Ow(III))- Experiment shows partial occupancy for each type 
of O w site; corresponding to an average of about 28 of the 
56 Ow sites shown occupied per large pore, (b) Way of occu¬ 
pying 28 Ow sites of (a) without any pair of water molecules 
approaching too close. Cu-Ow(IH) bonds shown in magenta, 
(c) Way of occupying 30 Ow-sites of (a). 



FIG. 3: Topology of water oxygens (Ow) in two different 
models for water structure of highly-hydrated Cu-BTC. (a) 
Model-28, corresponding to the Ow network in Fig. 2(b). 
(b) Model-42, a fullerene-like arrangement of 42 Ow- 


of about 28 H 2 0 molecules per large pore in arrange¬ 
ments that vary between pores. Physically, the partial 
occupancies are mandated because many of the potential 
sites are too close together to be simultaneously occu¬ 
pied. In Fig. §*>) and Fig. [2jc), we show two differ¬ 
ent choices for placing Ow that avoid unphysically close 
pairs of water molecules. The first has an arrangement 
of 28 water molecules with the Ow forming a polyhe- 
dra, and the second consists of 6 independent cyclic clus¬ 
ters of 5 water molecules, for 30 total water molecules. 
We refer to models with these starting configurations of 
Ow as “Model-28” and “Model-30”, respectively, where 
the number refers to how many water molecules per prim¬ 
itive Cu-BTC cell, or, equivalently, per large pore. 

The choice in Fig. [2jc) has the same topology as the 
C 2 s fullerene of T^ symmetry. [38] We show the topol¬ 
ogy in Fig. [3ja) . In fullerenes, each carbon has three 
carbon neighbors. In the equivalent water structure, 
a highly hydrogen-bonded structure is produced in an 
analogous way to the ice rules for three-dimensional ice: 
each Ow is covalently bonded to two hydrogens, and 
each Ow-Ow neighbor link has one hydrogen atom that 
is bonded to one Ow and forms a hydrogen bond to 
the other. These observations inspired us to look for 
other fullerene-like arrangements of Ow molecules that 


TABLE II: Comparison of symmetrized DFT and 
experimental^!] structures of dry Cu-BTC. Space group 
Fm3m with a = 26.2793 A. The standard deviation of the 
experimental structure refinement coordinates are 0.0001 or 
smaller. 




DFT 


Expt. 

Atom Wyckoff pos. 

x y z 

X 

y z 

Cu 

48(h) 

0.0000 0.2170 0.2170 

0.0000 

0.2164 0.2164 

O 

192(1) 

0.1833 0.2437 0.4478 

0.1831 

0.2438 0.4480 

C(l) 

96 (k) 

0.2035 0.2035 0.4308 

0.2037 0.2037 0.4302 

C(2) 

96 (k) 

0.1786 0.1786 0.3864 

0.1783 

0.1783 0.3865 

C(3) 

96 (k) 

0.1352 0.1352 0.2996 

0.1352 

0.1352 0.3000 

H 

96 (k) 

0.1189 0.1189 0.2654 

0.1202 

0.1202 0.2703 


might be accommodated in Cu-BTC. We looked at all 
low-energy fullerene geometries in the online Atlas of 
Fullerene Structures [38] from N = 20 to N = 60 at vari¬ 
ous orientations and scalings within the Cu-BTC large 
pore to see what structure would best accommodate 
preferred Ow-Ow distances of 2.9 A, while having one 
Ow near each exposed Cu. The best candidate struc¬ 
ture is a 42-molecule structure equivalent to the “C 4 2 
$45” fullerene, where the Ow have the topology shown 
in Fig. [3jb). We call this “Model-42”. For each model, 
as well as a model “Model-12”, with one H 2 0 per Cu, we 
initially randomized the hydrogen positions, then per¬ 
formed a full relaxation until all forces were converged 
within 0.03 eV A -1 . 


Results and Discussion 

The calculated structure of dry Cu-BTC is shown 
in Table [TT] and compared with a recent experimental 
X-ray powder diffraction refinement. |2lj Extremely good 
agreement (within 0.02 A) is obtained for all atomic po¬ 
sitions except for H, which are 0.14 A from the experi¬ 
mental positions. We note, however, that the experimen¬ 
tal refinement treated the organic ligand as a rigid body; 
and did not further refine the average hydrogen positions, 
which have little effect on the X-ray powder diffraction 
pattern in any case. 

As a bridge to the water cluster studies, we first investi¬ 
gated the binding energy of a single water molecule in Cu- 
BTC. We tested each Ow position found in the Wong-Ng 
al. (W-N) powder diffraction structure refinement. [2Tj 
For comparison, we also tested Ow hr the high-symmetry 
locations in Cu-BTC: the centers of the large, medium, 
and small pores, and the centers of the square and tri¬ 
angular windows between pores. The center of the small 
tetrahedral pore has been found to be a favored bond¬ 
ing site for certain sorbates in previous work on Cu- 
BTC. [2213S] We additionally investigated a medium pore 
“bonding”, site set 3.8 A from the center of the Cg 
ring, chosen because it is a favorable distance for water- 
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TABLE III: Calculated binding energies A E of a single water 3 . 
molecule in Cu-BTC. at various positions of the O w- “W- 
N” sites are from Ref. [2TJ Binding energies are broken down 
into van der Waals (vdW) and chemical (chem) contributions. 

Ow positions are given in crystallographic notation and all 
energies are in eV. 


O w position x y z A E vdW chem 

W-N I 0.1307 0.1307 0.1307 -0.14 -0.08 -0.06 

W-N II 0.2287 0.0364 0.0364 -0.12 -0.06 -0.06 

W-N III 0.1696 0.1390 0.0000 -0.47 -0.12 -0.35 

W-N III + full relaxation 0.1712 0.1413 0.0000 -0.53 -0.13 -0.40 U 

Large pore center 0.0000 0.0000 0.0000 -0.01 -0.01 -0.00 

Medium pore center 0.5000 0.5000 0.5000 -0.02 -0.01 -0.01 

Small pore center 0.2500 0.2500 0.2500 -0.12 -0.09 -0.03 

Square window 0.0000 0.0000 0.2564 -0.05 -0.03 -0.02 

Triangular window 0.1900 0.1900 0.1900 -0.18 -0.14 -0.04 

Medium pore “bonding” 0.4260 0.4260 0.4260 -0.06 -0.04 -0.01 






carbon separation. 140] For simplicity, the Cu-BTC frame¬ 
work was kept rigid in each case and only the hydro¬ 
gens were allowed to move, except in the case of W-N 
Ow (HI), where, in addition, a full relaxation was per¬ 
formed. We broke down the total binding energies A E 
into two parts: molecule-framework van der Waals en¬ 
ergy (vdW), and the excess part, due to chemical (bond¬ 
ing/antibonding and nondispersive electrostatic) interac¬ 
tions of the molecule with the framework (chem). The 
vdW (dispersion) energy was calculated using the ap¬ 
proximation of Grimme, only counting vdW interactions 
between the H 2 O molecule and the framework. The 
chemical component of the binding energy was obtained 
by simply subtracting the vdW contribution from the to¬ 
tal. The results are shown in Table ITTT1 

As expected water binds most strongly to the W-N III 
site (the one near the exposed Cu ion) and this bind¬ 
ing is mostly due to chemical interactions. The second- 
strongest binding site found is the triangular window, 
due to large van der Waals interactions, but this site is 
absent from the experimental structure refinements for 
some reason. The small pore center and W-N I and W-N 
II sites bind water molecules in the -0.12 eV to -0.14 eV 
range. Aside from the W-N III site, sites W-N I and W-N 
II have the strongest chemical interactions. Bonding is 
weak at the large and medium pore centers due to dis¬ 
tance from the framework and bonding is weaker in the 
medium pore than the large and small pores due mainly 
to the lack of chemical interactions with the framework. 

Now we give the results for Cu-BTC with multiple wa¬ 
ter molecules. The relaxed structures for the hydrated 
Cu-BTC models are shown in Fig. [4] Hydrogen bonds 
are shown as thin green lines. The key structural param¬ 
eters and energetics results are shown in Table |IV[ As 
the number of water molecules increases, the structure 
changes from individual water molecules bound to Cu to 
clusters of water molecules bound to one or more Cu, to 



FIG. 4: Structures of model water clusters inside large pore of 
Cu-BTC, relaxed via density functional theory, (a) Model-12; 
(b) Model-28; (c) Model-30; (d) Model-42, where the number 
refers to the number of water molecules inside the large pore. 
Hydrogen bonds are shown by thin green lines. 


TABLE IV: Key results for the different models for hy¬ 
drated Cu-BTC investigated with this work, compared with 
experimental data (where available). Distances D are in 
A and energy differences A E are in eV. Standard deviations 
of the experimental structure refinement results are given (in 
parentheses, with units of the least-significant digit). 


Study 

N w 

D Cu-Cu 

D cu-o 

Dcn-O ff 

<1 

AE/N W 

DFT (dry) 

0 

2.454 

1.932 




DFT Model-12 

12 

2.537 

1.946 

2.262 

-6.30 

-0.52 

DFT Model-28 

28 

2.582 

1.954 

2.241 

-16.66 

-0.59 

DFT Model-30 

30 

2.589 

1.950 

2.263 

-17.49 

-0.58 

DFT Model-42 

42 

2.586 

1.954 

2.217 

-26.81 

-0.64 

Expt. (dry) 

0 

2.498(1) 

1.930(1) 





Expt. (hydrated) 27.84 2.628(10) 1.951(6) 2.32(2) 
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FIG. 5: Structure of a 28-water cluster inside a medium pore 
of Cu-BTC, relaxed via density functional theory, (a) 100 
view (b) 110 view. 


a closed cage encompassing all of the water molecules. 
The trend in key interatomic distances with increasing 
numbers of water molecules agrees with experiment. 

Although Model-28 and Model-30 begin with atoms 
at positions suggested by experimental structure refine¬ 
ment, certain water oxygen positions relax as much as 
1.6 A, suggesting a significant (quasi)static contribution 
to the large experimental displacement factors found for 
water oxygens. [21] The relaxed DFT positions do remain 
within and near the inner surface of the large 13 A pore. 
Our Model-42 results show that 3.5 molecules per Cu are 
easily accommodated in this inner surface region. If wa¬ 
ter also occupies the interior of the large pore as well as 
the 11 A and 5.5 A pores, then Cu-BTC can accommo¬ 
date the 6.5 to 8.0 H 2 0 molecules per Cu ion suggested 
by experiment without requiring larger (defect) pores. 

The binding energy per water molecule for one H 2 0 at¬ 
tached to each Cu site is -0.52 eV, almost identical to that 
for an isolated H 2 0 molecule (Table III vs. Table IVI. 
In the highly-hydrated state configurations, the magni¬ 
tude of the binding energy per water molecule is even 
larger. Where does this enhanced stability come from? 
To determine this, we broke down the binding energy of 
each model into three parts: intracluster energy; cluster- 
framework van der Waals energy (vdW), and the excess 
part, due to chemical interaction of the cluster with the 
framework. The intracluster binding energy was calcu¬ 
lated by removing the framework and performing a DFT 
energy calculation on just the cluster. The remaining 
contribution to the binding energy was then broken down 
into vdW and chemical contributions as described above 
for the single molecule case. As a control, we also took 
a 28-molecule water cluster, moved it to the medium 11 
A pore, naming this model “Model-28MP”, fully relaxed 
the system (Fig. [5]), and broke down its energy into com¬ 
ponents. The results are shown in Table [Vj 

From Table [Vj the primary cause for the high stabil¬ 
ity of water clusters in Cu-BTC is the combination of 
intracluster (hydrogen bonding) interactions within each 
cluster and chemical interactions between the cluster and 
the framework, particularly when each exposed Cu site 
inside the large pore “sees” a molecule belonging to a 


cluster. The medium pores lack exposed Cu sites, and the 
Model 28-MP cluster in this pore has almost zero chem¬ 
ical interaction with the framework. Comparing Model 
28-MP and Model 28, the intracluster interactions are 
stronger in the medium pore, but the overall binding en¬ 
ergy of 28 water molecules is about 1.8 eV more stable in¬ 
side the large pore, Secondary influences on water cluster 
stability in Cu-BTC include van der Waals interactions 
and secondary chemical interactions. Secondary chemi¬ 
cal interactions results from two sources: enhancement of 
H 2 0-H 2 0 binding under an electric field |411 - fl5] (gener¬ 
ated by the exposed Cu ions), and additional H-bonding 
of H 2 0 H to O in the framework itself, seen in Fig. |4](b) 
and (d). 

All calculations were performed at zero temperature, 
and do not include thermal motion. Given that wa¬ 
ter is a liquid at room temperature, what can we say 
about the nature of the water in highly-hydrated Cu- 
BTC at room temperature? Here, we refer to a recent 
molecular dynamics simulation [44] on “MIL-53”, a Cr- 
based MOF with a smaller unit cell than Cu-BTC. They 
found, at room temperature, dynamic hydrogen bonds, 
one Ow (quasi)-statically bound to each exposed hy¬ 
droxyl unit, and additional water molecules that move 
fluidly. This is akin to the picture discussed in an NMR 
study on Cu-BTC,[2DJ where they conclude that there 
are two types of water molecules: one type bound to 
Cu, and the other fluid. The experimental X-ray powder 
diffraction refinement of Cu-BTC, showing 2.3 Ow bind¬ 
ing sites per Cu,[21] suggests that more water molecules 
are (quasi)-static than just those immediately bound to 
the Cu, but that any additional water molecules in the re¬ 
maining pore space may be fluid as their oxygen positions 
are not resolved. Clearly, further theoretical and experi¬ 
mental studies on the dynamics of water in Cu-BTC are 
needed. It would be useful to perform molecular dynam¬ 
ics studies of H 2 0 in Cu-BTC. The relevant time scale 
for hydrogen bonding rearrangements in H 2 0 is of order 
picoseconds (see e.g Ref. HH]) and measurements of liquid 
water diffusion ratesjlB] imply that Ow~0\v rearrange¬ 
ments occur on the scale of tens of picoseconds. These 
time scales are too long for ab initio molecular dynamics 
to be practical; classical molecular dynamics on a suit¬ 
able force-field model are required. On the experimental 
side, diffraction studies at low temperatures, where the 
water is completely frozen, would be particularly useful. 

Do our results provide guidance for the design of a 
MOF with weaker water binding, one that might be more 
useful for C0 2 capture from a moist gas mixture? In 
the following, we speculate that the distance separating 
the metal ions in a MOF with exposed metal ions is one 
factor that may affect the selectivity of H 2 0 versus other 
sorbates. For Cu-BTC hydrated with one water per Cu 
ion (Fig. ga)), the closest O^-Ow distance is 5.2 A, 
which coincides with the third-neighbor peak in the 0-0 
pair distribution function of water ice. m On the other 







6 


TABLE V: Breakdown of interactions of water clusters in Cu-BTC. Total cluster binding energy A E and cluster binding energy 
per water molecule (A E/Nw are broken down into three components: intracluster (intra), cluster-framework Van der Waals 
(vdW), and cluster-framework chemical interactions (chem). All energies in eV. 


Model 

N w 

A E 

total intra 

vdW chem 

A E/N w 
total 

intra vdW chem 

DFT Model-12 

12 

-6.30 -0.22 

-1.49 -4.59 

-0.52 

-0.02 -0.12 -0.38 

DFT Model-28 

28 

-16.66 -8.13 

-2.66 -5.86 

-0.59 

-0.29 -0.10 -0.21 

DFT Model-30 

30 

-17.49 -9.96 

-2.60 -4.93 

-0.58 

-0.33 -0.09 -0.16 

DFT Model-42 

42 

-26.81 -16.64 

-3.91 -6.25 

-0.64 

-0.40 -0.09 -0.15 

DFT Model-28MP 

28 

-14.81 -12.72 

-2.10 0.01 

-0.53 

-0.45 -0.07 0.00 


hand, an 0-0 distance of 5.8 A, is near a minimum of 
the 0-0 pair distribution function for both ice and liquid 
water. gum If one could synthesize a MOF with open 
metal sites such that the equilibrium Ow~Ow distance 
between water molecule bound to adjacent metal ions 
were 5.8 A, the formation of water clusters that bridge 
this unfavorable Ow~Ow distance might be inhibited, 
leading to a more favorable balance between CO 2 and 
water sorption. 


Conclusions 

Meta-GGA+U density functional theory calculations 
are shown to reproduce the experimental structure of 
dry Cu-BTC and also to accurately give the interaction 
energies of small water clusters. Applying these calcu¬ 
lations to models for hydrated Cu-BTC, we find that 
hydrogen-bonded water clusters have their stability en¬ 
hanced primarily by interactions with the exposed metal 
ions, and secondarily by van der Waals interactions, elec¬ 
tric field enhancement of water-water bonding, and hy¬ 
drogen bonding of water to framework oxygens. This ex¬ 
plains the great affinity of water for Cu-BTC and related 
metal-organic frameworks. 
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